Asymptotics of various partitions 
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1. Introduction. 

In [1] the author gave the first two terms in the asymptotic expansion of 
p(r, n) = the number of partitions of n into parts > r. In this paper we compute 
further terms until we get all digits of p(12, 1200). 

Further examples treated are "basic partitions" ( Nolan, Savage, Wilf [7]), 
3-colored partitions ( Kolitsch [5]), plane partitions strictly decreasing along 
rows ( Macdonald [6] ), the number of semisimple p-rings (Knopfmacher [4] ) 
and concave partitions (Snellman, Paulsen [8]). 



2. Partitions without small parts. 

We have the generating function 



l{x) = Y j p{r,n)x n - FU) 



where 



n=0 



oc 



fr-l(x) 



and 



n=0 j=l 

f r (x) = f[(l-x^ 
i=i 

In [1] we used derivatives of 

sinh(y / a£) 

to express the formulas. Here we will rather use associated Bessel functions 
with half integer index. By Example 5 in [1] we have 



fr(e- 1 ) = 



(I^IV-V, cxp (v B 23 B 2J+I (r) ) r _ B^t 

V^F P \f^ 2j(2j + l)\ j P V l ' 4 



Remember that formally we have 

,-a 

T(a) 
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Hence 



t a exp(rf- 1 + bt) = e bt ~t~ j+a 



J=0 J 



where 
Hence if 



OO n J. n j-j ^ 



£ = n + b 



then 



where 



$i(n) = i-^^c^. ^iUj 7_ (r+1/2+2j) ( 2v ^(2K) 

B 2 (r) 



For the second approximation belonging to the singular point x = — 1 of 
/ r (x) we obtain for even r 

fj-e-t) = ^-LL^-l)/2 eXP J ^ — + 22j=l 21(27+1)1 (^+1 (2) 



Assuming that the last part of the exponential function is 

00 

we find 

*2(n) = (-1)" ^^E^' (^J /-(,/2 + i/2 +2j )(v / «) 

where 

, S 2 (r) 
s = n > — 
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Higher order terms. 

We need to approximate f r (x) near x — cxp(2nih/k) where (h, k) = 1. This 
will lead to some formidable expressions containing Dedekind sums and various 
other sums of derivatives of cot(x) evaluated at x = jhn/k. We will save the 
reader from this. Instead we show how to do it numerically in an example. We 
continue Example 5 in [1]. Thus let r = 12 and n = 1200. Take k — 7 which is 
the first really difficult case. We consider the following three singular points: 

1. x = exp(27ri/7). 

In general we have 

[~kt fc(2)t _1 t 
F(exp(27ri/i/fc - *)) « J— exp tt a + ™(h, k) 
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where (h,k) — 1. Here 



S (h,k) = j2((jh/k))((j/k)) 

where 

{{x))=x-[x]- 1 -\tx$Z 
= otherwise 

Then 

,„ ^^ l~7t { <r(2) t 5iri 

since 

»d,7) = ^ 

Warning. 

In the beginning I used the table in Hardy-Ramanujan's famous paper [3] , 
but there are a lot of errors (even in the Chelsea edition of Ramanujan's collected 
papers ). So it is better to do your own computations. 



2. x — exp(4-7ri/7). 
Then 

J F 1 (exp(4 7 r J /7 - t)) w J ^- exp / 



2ir r 49i 24 14 



since 



3. x = exp(67ri/7) 
Then 




F(exp(27ra/7 - f)) « \ — exp , 

w 2tt 49i 24 14 



?(2) f iri 
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since 

*(3,7) = -l 

The other points with h = 4, 5, 6 are treated by conjugation. We have 

f r (x) = g{x)F{x) 

where 

n 

g(x) = - x 3 ) = 1 - x - x 2 + x 5 + x 7 - 2x 15 + ■ ■ -x 65 - x 66 



We expand 



numerically. Thus 



Put 



g(exp(2iri/7-t)) = J2 C ^ 

3=0 

c = 
ci = -56.249 - 246.445i 
c 2 = 3157.20+ 7835. 75i 

etc 



£ = 1200 - — 
s 24 



Then the contribution from x = exp(2ni/7) and x = exp(— 2-nijl) is 
10 . . o^ionnu r^o^ 3 / 4 +J'/ 2 



„ 7 A n f r 5ttz 2^12001 1 r ^(2) r /4+J/ % 

^ = 2 v ^ g * e cxp {u — — } } { w; 7 

= -26.49890 

Similarly we obtain 

$7,2 = 9.34180 
$ 7 ,3 = 8.86403 

and 

$ 7 = $ 7j i + $7,2 + $7,3 = -8.29306 
Taking 16 terms we get 

$i= 4 90015 90791 72948 38443 67842.13295 

$ 2 = 33 28709 71879.70797 
$ 3 = 124 91942.03287 
$ 4 = 51637.99481 



$5 = 


: 661.16304 


$ 6 = 


= 164.15289 


$7 = 


= -8.29306 


$ 8 


= 3.22670 


$g = 


= -0.19405 


$10 


= 0.23922 


$11 


= 0.00072 


$12 


= 0.34102 


$13 


= 0.89399 


$14 = 


= -0.64967 


$15 


= 0.94330 


$16 


= 0.44736 



Summing up 

16 

$j ■■= 4 90015 90791 72981 67278 84124.17000 

3=1 

which agrees with 

p(12, 1200) = 4 90015 90791 72981 67278 84124 



Remark. 

We have $n = 0.00072 but from k = 12 the $& become much larger. This 
depends on the fact that the constant term 

n 

IJ(l-exp(27ryty*)) 

3=1 

is no longer zero. 



3. Basic partitions. 

Our next example concerns a limiting case of" basic partitions" (see Nolan, Savage, Wilf 
[7] ). Thus we have the generating function 

00 1 _1_ 3 00 

/(*) = II ^ = 

7 = 1 n=0 
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Then 



and 



which gives 



where 



=1 ( X - xJ ) 2 F ^ 2 ) 



/(exp(-i)) 



F(exp(-t)) 2 _ 1 fl/2 
F(exp(-2i)) ~ 2^H 

1 /a\3/ 4 



t^exp(-) 



General k. 
Only a; 



exp(2Trih/k) with odd fc are singular points of /. We have 

F(cxp(2Trih/k-t)) 2 



f(exp{2irih/k - i)) 
Vkl 



Put 



and 



Then 



where 

A(l,n) = 
A(3,n) = 
A(5,n) = 
A(7, n) = 
A(9,n) = 

1(11, n) 
1(13, n) 



F(exp(47ri/i/fc - 2t)) 
- 2 ^ exp {iri(2s(h, k) - s(2ft, fc))} exp {p^j 

T(ft,fc) = 2s(h,k) - s(2h,k) 

A(k, n) = ^ CX P {niT(h, k) - 2mnh/k) 
(fc,fe)=i 

t / \ 1 /fc" / a \ 3 / 4 -., , T ,„ /an. 



1 

2 cos(27rn/3 — 7r/6) 

2 {cos(2t™/5 - 2tt/5) + cos(4t™/5 - tt/5)} 

2 {cos(27rn/7 - 9tt/14) + cos(4?m/7 - tt/14) + cos(6tto/7 - 3tt/14)} 
2 {cos(2t™/9 - 8tt/9) + cos(4t™/9 - 4tt/9) + cos(87rn/9 - 2tt/9)} 

f cos(27rn/ll - 25tt/22) + cos(47m/ll - 7tt/22) + cos(6vrn/ll - tt/22) 
~ \ + cos(87rn/ll - 9tt/22) + cos(10tt/11 - 5tt/22) 

cos(27rn/13 - 18tt/13) + cos(4ttti/13 - 9tt/13) + 
2 { cos(6t™/13 - 6tt/13) + cos(87rn/13 + 2tt/13)+ 
cos(10t™/13 - tt/13) + cos(127m/13 - 3tt/13) 



Example. Let n = 1000. Then 

amoo = 17293 58313 74933 37582 44155 69812 30246 17584 
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We have 

$1 = 17293 58213 74933 37582 44155 69812 30245 67360.51591 
$ 3 =0 

$ 5 = 50530.87514 
$ 7 = -305.16980 

$9=0 

$n = -1.78371 
$13 = -0.43964 



$i + • • -$i3 = 17293 58213 74933 37582 44155 69812 30246 17583.99790 
so the error is —0.0021. 

4. 3-colored partitions. 

We have the generating function 

(1 - x^f 



/Lr: H:l ,^;1 



n=0 



This case caused me a lot of trouble until I realized that one has to use the 
exact functional equation for 

F( X ) = f[(i- x n- 1 

3=1 

namely 

T f 7T 2 t 1 „, f 4rr 2 



nexpH)) = ^-exp|---|F(exp|- — 

Now 

_ F(xfF(x3) 

/W- F(x 9)3 

gives 



/(exp( - t} } = V 4^ exp { Yt + 1 } { 1 - 3 cxp { - %■ 

The other terms coming from i^(exp(— 47r 2 /t)) and -F(exp(— 47r 2 /3i)) give terms 
of the form exp(— c/t) where c > which transform into ordinary Bessel func- 
tions J_3/2(2-y/cJ) which are very small for large £. 
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We get as usual with £ = n + 7/8 
More generally we have 

F(exp{2irih/k - t)) = 

- exp ^zs(h, k) + m - -| F(exp(|2**& /* - j) 
where 

/ift/ = — 1 (mod fc) 

In particular for k = 2 we get 



243tt ^ 8i 8 72i 8 



and thus 

-<»> - <-«"7asr { (sf '- <2 ^?> + 3 (wf '-^'Vf >' 

Further computations give (one has to be careful when k is divisible by 3 ). 
$ 3 (n =2cos )-== — /_3/ 2 (2W— ) 

.2 \ 3/4 



$ 4 (n =2cos h-)\/ J-3/2(2\/— ) 

1 ; V 4 8 7 V 484tt ^32^ 3/21 V 32 ; 

= 2 | C0S( — - T> + C ° S( — - 5 } | V 486^ \M J 7 - 3/2(2 V MT } 
, . „ ,27m 5tt, 1 / tt 2 \ 3/4 r . 

_ . . „ f ,27rn 117T. ,47m tt . ,67m 137T. 1 

* r (n) = 2 |cos(— - — ) + cos(— - -) + cos(— - — ) j 



486tt V98^J 7 3/2 V 98 ^ 

Numerical example. 
Take n = 200. Then 

a 200 = 1747 47949 05123 77771 22300 

and 
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$! = 1747 47949 05123 51995 85303.8855 

$ 2 = 25778 34263.0509 

$ 3 = -2 97748.8532 

$ 4 = 474.3640 

$5-0 

$ 6 = 5.5359 

$ 7 = 0.5261 

Thus 

$i + • • $ 7 = 1747 47949 05123 77771 22298.5093 
so the error is —1.4907 



5. Plane partitions strictly decreasing along rows. 

We have the generating function (see Macdonald [6] p. 83) 

OO OO 

f( X ) = U(1- xiy^W =J2"nX n 
j = l n=0 

Consider 

OO 

g(t) = log(/(cxp(-i))) = - J2 3 {Iog(l - exp(-(2j - l)i)) + log(l - exp(-2ji))} 

3=1 

OO OO 

= J2 J2 i {exp(-fc(2j - l)t) + exp(-2jfci)} 

j = l fe=l 

We form the Mellin transformation of g(t) 

/• OO 

g{s)= / g(t)r s dt = 
Jo 



OO OO ^ 

rooEE{ 

i=lfc=l k 



k 1 + 8 (2j - l) s /c 1 + s (2j) s 



Ir( s K(i + s){^-i) + (i-2- s K( s )} 



5 (is 



To compute the inverse Mellin transformation 

Y yc+ioo 

9{t) = 2^ / ?(s)r 
(c > 2) we need the residues 

iMs(*rv = 2) = ^ 
i? es (g( S )t- s , s = i) = ^ 
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^ Wt -,- = 0) = I{^(-l)-^ + ^ 

Re S (g( S )t- s ,s = -l) = ^ 

Res(g(s)t- S ,s = 1 - 2k) = for k > 1 

1 t 2fc 
ResQ( S )t- s , s = -2k) = - 2fc)?(-l - 2fc)^y for fc > 1 



Hence 

/(exp(-i)) 



Use the notation 

2 

24 



c = 2- 1 / 4 cxp(^(-l)/2) 



Remark. 

?'(— 1) is known as Kinkelins constant (see [2] for some history ) 
First we replace 



I 1 ^ ?(1 - 2fck(-l - 2fc)i 2fe i _ ^ 

° p l€ <»>■ }= 1+ £ c " 

by 1. We also disregard exp(i/48) which is merely the translation n — > ra+1/48 : 

e 

/(exp(-t)) w ci 1/24 exp(a^ 2 + foi^ 1 ) = 



oo 

fc 



°° °° a J'5fei-2j"-fe+l/24 ^ ~ a j b kt2j+k-l/24-l 

C ' 



f-; ^ " c j!fc!r(2j + jfc - 1/24) 

3=0 k=0 J 3=0 k=0 ■' W 7 ; 



Hence 



(, \ 25/48 oo -. , C\i 
l) E^(y) ^M 2 ^) 

We also have to consider terms coming from c 2 t 2 + + • • • ■ We have 

1 



C2 



c 4 



5760 
313 



464486400 
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91207 
° 6 ~ ~ 8026324992000 

etc 

The term coming from c 2 t 2 will be 

-, / i \ 25/48+1 oo -, / e\3 

Let us take n = 200. Then 

a 200 = 233 1 86 5 1 62 1 79 6 3 5 3 6 5 70 1 4 
Taking eight terms we get 

$1 = $1,0 + $1,2 H $1,14 = 

233 18651 62179 61259 31636.4295 

so we get 14 correct digits out of 23. There is a lot of cancellation so we have 
to compute with many more digits than 23. 



Second term. 
Let 

<&(*)= log(/(-exp(-t)) = 

OO 

3 {log(l - exp(-(2j - l)t)) - log(l - exp(-2(2j - l)t)) - log(l - exp(-2ji))} 

with Mellin transformation 

g 2 {s) = \t{ S ),{1 + a) {(5 • 2" s - 2 • 2~ 2s - - 1) - (1 - 2- s )^( S )} 
with residues 

Re S (g 2 (s)t- s ,s = 2) = ^ 
Res{g 2 { S )t- s ,s = l) = -^ t 
Res( M s)t-°,s = 0) = <> { -l )+ 1 ^ + 1 ^ 

Res(g 2 (s)t- S ,s = -1) = ^ 

13f 2 

Res(g 2 (s)t- s , S = -2)- ' 



i?e S (g 2 (s)t- s , S = -4) 



5760 
433t 4 
1451520 
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Res(g 2 {s)t~ s , s = -6) 



97091200 

etc 



Hence 



exp 

We get 



13^ 433i 4 7873r 



5760 1451520 87091200 



/, \ 13/24 oo . , j _ 

$ 2 , (n) = (-l)«2 1 ^ exp( ,' ( _ 1)) M £ J_ |i I-^-is/iapVk? 



where 

_C(3) 



6 2 = - — 
2 48 

£ = n H 

^ 48 

We obtain & 2 ,k for fc = 2, 4, 6, ■ ■ ■ by expanding 

13t 2 



exp 

The result is 



5760 

$ 2 = 2281 49810.1011 



and 

$i + $ 2 = 233 18651 62179 63540 81446.5432 
so we get 17 correct digits out of 23. 



6. Semisimple p-rings. 

Let s(n) = the number of isomorphy classes of semisimple p-rings of order 
P n (P =prime, see Knopfmachcr [4], p. 64). We have the generating function 

oc oo oo 

n=0 j=l k=l 

To get the first approximation we consider 

OC OC OC / - i 9 , \ 

g(t) = log(/(exp(-t))) = E E E 

j=l k=l r=l 



12 



The Mellin transformation is 

g(s)= / gity^dt = r(s)?(l + s)?(s)?(2a) 



with residues 

Res(g(s)t- S ,s= 1/2) = 1^(1/2)^(3/2) 



Hence 



Let 



Then 



9(() = «_^ + £ (2) ! + ^ (1/2W3/2rl/2 



a = ^(2) 2 
&=^?(l/2M3/2) 

/(exp(-t)) = (2 7 r) 3 / 4 t- 1 /4 cxp ( ai -i + w -i/2j 

2 j])k t -j-k/2-l/A 



oo oo 



(2-) 3/4 EE 



j=0 fc=0 J 



oo oo 



^^ j!fc!r(j + fc/2 + l/4) 



Hence 



oia /a\3/8^, 1 / /n\V 4 \ 
$i(n) = (2tt) 3 / 4 (-J Em Ha) ^-^xM 

To get the second term we consider 

oo oo 

g 2 (t) = log(/(exp(-t))) = - E E - (- 1 ) J ' £2 CX P(-^ 2 *)) 

j=0 fc=0 

After some computations one obtains 

g 2 (s) = r(a)?(l + «)?(*) {-1 + 3 • 2- s + 2- 2s - 3 • 2~ 3s + 2~ 4s } 
with residues 
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Res(g 2 ( S )t- s , s = 1/2) = {3^ - l} ?(l/2) ? (3/2) = -| 

^ /- / n , ^ 31og(27r) log(t) 
Res(g 2 ( S )t- s , 8 = 0) = - 

so 

/(-exp(-t)) - (27r) 3 / 4 t- 1 / 4 exp{a 2 i- 1 +M" 1/2 } 
which is exactly /(exp(— t)) with a — > a 2 and b — > 6 2 - Hence 

3/8 00 / / \ 1 / 4 \ fe 

* 2 (n) = (-l)«( 27 r) 3 / 4 (£) £ 4/2-3/4(2^) 



j=0 



Numerical Example. 
Let n = 200. Then 



s(200) = 26122 95856 86401 



while 

$i(200) = 26122 95856 68813.6838 
$ 2 (200) = 290 18442.4542 

$1 + $ 2 = 26122 95852 86711.2679 

so we get 9 correct digits out of 15. The largest terms in $1 are of size 10 22 so 
there is a lot of cancellation. Knopfmacher gives the estimate 

7T 2 

log(s(n)) ~ — 

We can improve this by using the asymptotic formula for Bessel functions. We 
have for large z 

Ir{z ) „ CX P( Z ) L _ Ar2 - 1 + ( 4r2 - !)( 4r2 - 9 ) 



f%Tz I Sz 2\(Sz) 

Putting 

z = 

00 

_ k 3 

T ~ 2 ~ 4 



and summing over k we obtain 



. , (2ttV^) 1/4 L ^- , /n\V4 6 2 1 
s(n) ~ v /l , exp <^ 2^/on + b - - — ^ 
V ' V2n 5 /8 p \ v W 16a 



For n = 200 the right hand side is 0.28 • 10 15 so we almost get one correct digit. 
7. Concave partitions 
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We have the generating function ( see [8] ) 

oc oo 

/(x)=n(i-^° +i)/2 )- i =E a ^" 

j=l n=0 

Observe that a n is also the number of partitions of n into triangular numbers. 
Then 

g(t) = log(/(exp(-f)) = £ £ I exp(-M^±l)!) 
j=i k=i 

and the Mellin transformation 

1 



J= ltlfc 1+s (j(.7 + l)) s 
= 2 s T(s)C(l + s)H(s) 

where 

We have to expand H(s) near its poles and the poles of r(s)£(l + s). Now 



H(s) = «2s) + f: -tja + j)- 8 - 1 } 



fe=l ^ ' 

We want the residues at the singular points of 



k=i 



G( s )-2 s r( s )c(i + 5 r s |c(2.s) + E( k )C(2 5 + fc) 

The result is 



since 



Res(G(s),s = 0) = log(^) - 7 + £ = lo §(^) 

3=2 3 



\J : = 7 = Euler's constant 

Res(G(s),s = --k) = 16 J fc(fc + 1) , 
15 



We have 

Res(G{s),s = -2k) = 

which is by no means trivial. It follows from the following identity: 
Proposition. (Joakim Petersson). 

Hence we have 



/(exp(-t)) = ^ exp _ f (-) • rV» + £ 8^(fc+T)! 



fe=0 

Put _ 

V^F 3 

The main part of the first approximation of a n comes from 

47r exp(«t J-4r^ > tel. m ft-)) 



9o 



(here we have l/r(0) = 0). To get a more precise estimate we expand 



and get 

_ c k ~ gini/ 2 - 2 -"/ 2 

qk 47r ^ - 1 - fc / 2 ) 

We take a numerical example: Let n = 2000. Then 

02000 = 27 79195 55391 39291 



and 



Thus 



qo = 28 36990 09214 40082.213 
qi = -58382 21952 56521.094 
q 2 = 594 72145 84992.877 
q 3 = -8 07564 32951.404 
<74 = 10135 44447.136 

$1 « 9o + Qi + 92 + 93 + 94 = 
27 79194 61978 60047.307 
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so we get almost 7 correct digits. 
Remark. 

In [8] they show that 

log(a„) ~ An 1 / 3 

as n — > oo. Using the saddle point method and our /(exp(— t)) we obtain 

A=^r^ V/3 C( 3 )2 /3 
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